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Abstract 

We derive an analytical expression for the scattering of a scalar wave from 
a perfectly conducting self-affine one dimensional surface in the framework 
of the Kirchhoff approximation. We show that most of the results can be 
recovered via a scaling analysis. We identify the typical slope taken over 
one wavelength as the relevant parameter controlling the scattering process. 
We compare our predictions with direct numerical simulations performed on 
surfaces of varying roughness parameters and confirm the broad range of 
applicability of our description up to very large roughness. Finally we check 
that a non zero electrical resistivity provided small does not invalidate our 
results. 

I. INTRODUCTION 

Although studied for more than fifty years [1] wave scattering from rough surfaces re- 
mains a very active field. This constant interest comes obviously from the broad variety of 
its applications domains which include remote sensing, radar technology, long range radio- 
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astronomy, surface physics, etc., but from the fundamental point of view, the subject has 
also shown a great vitality in recent years. One may particularly cite the backscattering 
phenomena originating either from direct multiple scattering [2,3] or mediated by surface 
plasmon polaritons [5-9]. Remaining in the context of single scattering a large amount 
of works have also been devoted to the development of reliable analytical approximations 
[10-13]. In all cases, the efficency of any analyitical approximation relies on a proper de- 
scription of the surface roughness. In most models the height statistics are assumed to be 
Gaussian correlated. In this paper we address the question of wave scattering from rough 
self-affine metallic surfaces. Since the publication of the book by B. B. Mandelbrot about the 
fractal geometry of nature [14], scale invariance has become a classical tool in the description 
of physical objects. In the more restricted context of rough surfaces, scale invariance takes 
the form of self-affinity Classical examples of rough surfaces obeying this type of symmetry 
are surfaces obtained by fracture [15] or deposition [16]. More recently it was shown that 
cold rolled aluminum surfaces [17] could also be successfully described by this formalism. 
When dealing with wave scattering from rough surfaces, this scale invariance has one major 
consequence of interest, it is responsible for long range correlations. After early works by 
Berry [18], lots of works have been performed to study the effects of fractal surfaces on 
wave scattering. Most of these works were numerical (see for example Refs. [19-27]) and 
very few analytical or experimental results have been published. Notable exceptions are 
due to Jakeman and his collaborators [28,29] who worked on diffraction through self-affine 
phase screens in the eighties and more recent works applied to the characterization of growth 
surfaces [30-32]. We recently gave a complete analytical solution to the problem of wave 
scattering from a perfectly conducting self-affine surface [33] in the Kirchhoff approxima- 
tion. In the following we present a complete derivation of this expression and we deduce 
from it analytical expressions for the width of the specular peak and the diffuse tail. These 
results are compared to direct numerical simulations. We show evidence that the crucial 
quantitative parameter is the slope of the surface taken over one wavelength. 
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II. THE SCATTERING SYSTEM 



The scattering system considered in the present work is depicted in Fig. 1. It consists 
of vacuum in the region z > ((x) and a perfect conductor in the region z < ({x). The 
incident plane is assumed to be the xz-plane. This system is illuminated from the vacuum 
side by an s-polarized plane wave of frequency u = 2ir/\. The angles of incidence and 
scattering respectively are denoted by 9 and 9, and they are defined positive according to 
the convention indicated in Fig. 1. 

In this paper we will be concerned with 1 + 1-dimensional self-afline surfaces z = ({x). 
A surface is said to be self-afline between the scales £_ and £+, if it remains (either exactly 
or statistically) invariant in this region under transformations of the form: 



Ax — > fiAx, 



(la) 



(lb) 
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FIGURES 




FIG. 1. The scattering geometry considered in this paper. 

for all positive real numbers /i. Here H is the roughness exponent, also known as the 
Hurst exponent, and it characterizes this invariance. This exponent is usually found in 
the range from zero to one. A statistical translation of the previous statement is that the 
probability p(A(; Ax) of having a height difference in the range [A(, A( + dA(] over the 
(lateral) distance Ax is such that: 



p (AC; Ax) dA( = p U H A(; ptAx) dpi" A( 



H 



(2) 



Simple algebra based on the scaling relation (1) gives that the standard deviation of the 
height differences ((x + Ax) — ((x) measured over a window of size Ax can be written as 



a(Ax) 

and the (mean) slope of the surface as 

s(Ax) = 



9l-H xH 



Ax n , 



(3a) 



l-H 



Ax 



(3b) 



In these equations I denotes a length scale known as the topothesy. It is defined as a {I) = I 
(or s(£) = 1). 

Alternatively, Eq. (3a) can be written in the form 



a(Ax)=a(X) (^-) H = Xs(X) 



A 



A 



(4) 



where we use the wavelength A of the scattering problem as a normalization length. Here 
cx(A) and s(A) are respectively the typical height difference and slope over one wavelength 
as defined by Eqs. (3). Note that we could have used any length scale for the normalization, 
like for instance the topothesy. However, the choice made here was dictated by the physical 
problem studied. Using similar scaling arguments one can show that the power density 
function of the height profile V(k) depends on the wave number A; as a power law: 



V(k) 



/oo 
((x) exp(ikx)dx 
-oo 



2 oc k' x ~ 2H . (5) 



In the case of a Gaussian height distribution, the probability p(A(; Ax) reads: 

The self-affine profile is thus fully characterized by the roughness exponent H, the slope s(A) 
(which is nothing but an amplitude parameter) and the bounds of the self-affine regime £_ 
and £+. 

Numerous methods have been developed to estimate these parameters (see for example 
Ref. [34]), most of them use the expected power law variation of a roughness estimator 
computed over spatial ranges of varying size. This roughness estimator can be a height 
standard deviation, the difference between the maximum and the minimum height, etc. It 
is also classical to use directly the power density function of the profile. More recently the 
wavelet analysis has beeen shown to offer a very efficient method to compute the roughness 
exponent of self-affine surfaces [35]. 



III. SCATTERING THEORY 

In the following we consider the scattering of s-polarized electromagnetic waves from a 
one-dimensional, random, Gaussian self-affine surface z — ((%)■ It will be assumed that the 
lower limit of the self-affine regime £_ is smaller than the wavelength, A, of the incident wave. 
For the present scattering system, where the roughness is one-dimensional, the complexity 
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of the problem is reduced significantly. The reason being that there is no depolarization 
and therefore the original three-dimensional vector scattering problem reduces to a two- 
dimensional scalar problem for the single non- vanishing 2nd component for the electric field, 
&(x,z\iv) = E y (x,z\u>), which should satisfy the (scalar) Helmholtz equation 

with vanishing boundary condition on the randomly rough surface z = ((x), and outgo- 
ing wave condition at infinity. In the far field region, above the surface, the field can be 
represented as the sum of an incident wave and scattered waves: 

®(x,z\u) = %(x,z\uj) + [°° p-R( q \k) e iqx+iaoMz , (8) 

J-oo 2n 

where the plane incident wave is given by: 

$o(x, z\uS) = exp {ikx — ia (k, co)z} (9) 
and R(q\k) is the scattering amplitude. In the above expressions, we have defined 



u- 2 



a (g,^) = y(^-J -<? 2 , $la (q,uj) > 0,%a (q,uj) > 0. (10) 

Furthermore, the (longitudinal) momentum variables q and k are in the radiative region 
related to respectively the scattering and incident angle by 

k = — sin 9q, (11a) 

c 

g = -sin#, (lib) 

c 

so that the ^-components of the incident and scattering wavenumbers become 

a {k, uj) = — cos#o, (H c ) 

a (q, uj) = — cos#. (Hd) 

c 

The mean differential reflection coefficient (DRC), also known as the mean scattering 
cross section, is an experimentally accessible quantity. It is defined as the fraction of the 
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total, time-averaged, incident energy flux scattered into the angular interval (9,9 + d9). It 
can be shown to be related to the scattering amplitude by the following expression [36]: 



<9/L \ 1 tu cos 2 1 



l^(#)| 2 ). (12) 



\ 89 I L 2nc cos^o 
Here L denotes the length covered by the self-affine profile as measured along the x-direction, 
and the other quantities have been defined earlier. The angle brackets denote an average over 
an ensemble of realizations of the rough surface profiles. Moreover, the momentum variables 
appearing in Eq. (12) are understood to be related to the angles 9 and 9 according to 
Eqs. (11). 

We now impose the Kirchhoff approximation which consists of locally replacing the sur- 
face by its tangential plane at each point, and thereafter using the (local) Fresnel reflection 
coefficient for the local angle of incidence to obtain the scattered field. Notice that dealing 
with a surface whose scaling invariance range is bounded by a lower cut-off £_ does ensure 
that the tangential plane is well defined at every point. Within the Kirchhoff approximation 
the scattering amplitude can be expressed as [36]: 

R( q \k) = ^— f L/2 dx e- iqx - iao{q ^ {x) Mv(x\uo), (13a) 

2a (q, uj) J-l/2 

where Afo(x\u) is a source function defined by 

jV (x\u) = 2d n $ (x,z\u)\ z=ax) . (13b) 

Here d n denotes the (unnormalized) normal derivative defined as d n = —('(x)d x + d z . 

By substituting the expression for the scattering amplitude , Eq. (13a), into Eq. (12), one 
obtains an expression for the mean differential reflection coefficient in terms of the source 
function J\f Q (x\u); the normal derivative of the total field evaluated on the rough surface. 
After some straightforward algebra where one takes advantage of the fact that the self-affine 
surface profile function ((x) has stationary increments, one obtains the following form for 
the mean differential reflection coefficient 

/dR s \ uj 1 fcos[(9 + 9 )/2}\ 2 f L / 2 



89 I 2-nc cos # \cos [(9 - 9 ) /2] J J-l/2 



J dv exp ji— (sin# — sin# )f j Vt(v ), (14a) 



where 

n(y) = (exp {-i^[cos0 + cos# ] AC(f)|^ , (14b) 

with A((v) = ((x) — ((x + v). Note that the statistical properties of the profile function, 
((x), enters Eqs. (14) only through Q(v). With the height distribution p(5(; Sx) introduced 
earlier, Eq. (6), one may now analytically calculate the ensemble average contained in Q(v). 
For a Gaussian self-affine surface one gets 

Q(v) — J dz exp j— i— (cos 6 + cos# ) p( z ] v ) 

( uj cos 9 + cos 6>n , , . , i it tt\ , 
!X P <! ( - ^ -s(X)X v J }. (15) 



By in Eq. (14) making the change of variable 



u = v 



uj cos 6 + cos #o 



s(A)A 



(16) 



and letting the length of the profile extend to infinity, L — > oo, one finally obtains the 
following expression for the mean differential reflection coefficient: 
ldR s \ s(A)~^a-(*- 1 ) cos^ /^tan 



90/ v^costfo cos 3 \as- _1 s(A)5- 



^ -V^TTfV » ( 17a ) 



where 



a = 27rv / 2cos -7- cos — — -, (17b) 



oo 



£ a ( x ) = -*-/" rfA;e fa e- |fcr . (17c) 



oo 



and (0 < a < 2) 



2tt 

The quantity £ a (x) is known as the centered symmetric Levy stable distribution of index 
(or order) a [37]. This distribution can only be expressed in closed form for some partic- 
ular values of a —; a — 1 and a = 2correspond to the Cauchy-Lorentzian and Gaussian 
distributions respectively, C1/2 and £1/3 can be expressed from special functions. When the 
a-index in the Levy distribution C a (x) is lowered from its upper value a = 2 (Gaussian 
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distribution), the resulting distribution develops a sharper peak at x = while at the same 
time its tails become fatter. It is interesting to note from Eqs. (17) that the wavelength, 
A = 2-kc/uj, only comes into play through the slope s(A). The behavior of the scattered 
intensity is thus entirely determined by this typical slope s(A) and the roughness exponent 
H. 
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FIG. 2. The mean differential reflection coefficient, (dR s /dd), vs. scattering angle, 6, for a 
perfectly conducting self-affme surface. The plotted curves are the prediction of Eqs. (17). The 
Hurst exponent in all cases are H = 0.7 and topothesies I are ranging from I = 10~ 2 A (S(X = 0.016) 
down to 10~ 6 A (s(A) = 0.25) as indicated in the figures. The angles of incidence were (a) 6>n = 0° 
and (b) 9 = 0°. 

In Figs. 2 we show the mean differential reflection coefficient as obtained from Eqs. (17) 
for Hurst exponent H = 0.7 and different values of the slope s(A) ranging from 0.016 to 
0.25. The angles of incidence were 9 = 0° (Fig. 2a) and 50° (Fig. 2b). It is observed from 
these figures that as the amplitude parameter s(A) is decreased, while keeping the other 
parameters fixed, the portion of the scattered intensity scattered diffusely is reduced, while 
the power-law behavior found for the non-specular directions survives independently (within 
single scattering) of the amount of light scattered specularly. Furthermore, as the Hurst 
exponent is decreased (results not shown), and thereby making the topography rougher at 
small scale, the mean DRC gets a larger contribution from diffusely scattered light. This is 
a direct consequence of the properties of the Levy distribution mentioned above. 
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In order to discuss the features of the mean DRC which can be seen in Figs. 2 we will 
now proceed by discussing the behavior of the specular and diffuse contribution to (dR s /d0), 
i. e. close and far away from the scattering angle 9 = 9q respectively. 



A. The specular contribution 



We start by considering the specular contribution to the mean differential reflection 
coefficient. This is done by taking advantage of the asymptotic expansion of the Levy 
distribution around zero [38] 



c a (x) = — r (- 

ircx \a 



2r U) 



X 



+ 0(x A ). 



(18) 



By substituting this expression into Eqs. (17) one finds that the mean DRC around the 
speular direction 9 = 9 should behave as follows (59 < 1) 



'OR, 
09 



1 

2H 



x 



1 + w lz^ tanSo+ W 



2H 



4 



V 



w tan 9 



1 \ 2H> 



T{±j) (2«s/27rcos0 o ) 

From this expression it follows that the amplitude of the specular peak should scale as 
ldR s ' 



89 



^ 2H ' 



0=0o 2V2~nH(2V2n cos 9 ) H 1 s(A) 1 /^' 



(20) 



and that the peaks half width at half maximum, w, should be given by 
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FIG. 3. The full single scattering solution (solid line), Eqs. (17), for the mean differential 
reflection coefficient vs. scattering angle for a perfectly conducting self-affine surface compared 
to its specular (dotted line) and diffuse expansions (dashed line) as given by respectively Eqs. (19) 
and (24). The surface parameters used were H = 0.7 and £ = 10 _4 A (s(X) = 0.063), and the light 
was incident normally onto the rough surface. 



w(H,s(X),9 ) 



\ WzflT (2^71- cos # ) "~ l s(X) l/H . (21) 



It is worth noting that in the above expression the width of the specular peak depends 
on the wavelength A via the typical slope over one wavelength s(A). In case of Gaussian 
correlations, there would have been no dependence on the wavelength, the peak width w 
being simply proportional to the ratio cr/r, RMS roughness over correlation length. 

In order to test the quality of the specular expansion, Eq. (19), we show in Fig. 3 a 
comparison of this expression with the full single scattering solution obtained from Eqs. (17) 
for a surface of roughness exponent H = 0.7 and of slope over the wavelength s(A) = 0.063 
(£ = 10 _4 A) in case of normal incidence. The amplitude of the specular peak is seen to be 
nicely reproduced, but this expansion is only valid within a rather small angular interval 
around the specular direction 9 = 9 . 

It is interesting to notice that in the case of a non-zero angle of incidence, 9q 7^ 0°, the 
specular peak is slightly shifted away from its expected position 9 = 9 due to the presence 
of a non-vanishing term in Eq. (19) linear in 59. In this case the apparent specular peak is 
located at 9 = 9 + A9 , where A9o (A9 ~ w 2 < id) scales as 

Aflo ^ X tanfl (2y/2n cos 9 ) ^ s(\) 2 / H = ^L-± ta n9 w 2 (H, s(\),9 ) (22) 

H 1 (^) ^ 1 AH 

Such a shift has not, to our knowledge, been reported earlier for non self-affine (or non 
fractal) surfaces. Hence, due to the self- affinity of the random surface, we predict a shift, 
A# , in the specular direction as compared to its expected position at 9 — 9 . Notice that 
this shift vanishes for a Brownian random surface (H = 1/2). Moreover, for a persistent 
surfaces profile function (H > 1/2) the shift is positive while it becomes negative for anti- 

11 



persistent profile (H < 1/2). Unfortunately the specular shift A9 is probably too small to 
be observable experimentally for realizable self-amne parameters. 



B. The diffuse component 

We now focus on the diffuse component to the mean differential reflection coefficient, 
i.e. the region where 9 is far away from 9 = 9 . Now, using the expansion of the Levy 
distribution at infinity (the Wintner development) 



„ . , r(l + a) . fan\ / 1 \ 



we get the following expression for the diffuse component of the mean DRC {9 ^ 9 ) 

3-2H 

W • (24) 



dR s \ _ T(l + 2H) sm(irH) s(A) : 



cos^a 



09 I (Ait) 1 "- 1 cos9 



sin 



2 

In Fig. 3 the above expression is compared to the prediction of Eqs. (17). We observe an 
excellent agreement for angular distances larger than ten degrees. Moreover, it should be 
noticed from Eq. (24) that the mean DRC is predicted to decay as a power-law of exponent 
— 1 — 2H as we move away from the specular direction. For smooth surfaces (corresponding 
to small values of s(A)) this behavior results directly from a perturbation approach where 
the scattered intensity derived directly from the power density function of the surface. As 
shown above, in the case of self-affine surfaces the latter is a power law of exponent — 1 — 2H. 
Our results extend then the validity of this power law regime to steeper surfaces. 



IV. SCALING ANALYSIS 

It is interesting that most of the non-trivial scaling results derived above can be retrieved 
via simple dimensional arguments. Let us examine the intensity scattered in direction 9; 
in a naive Huyghens framework two different effects will compete to destroy the coherence 
of two source points on the surface i) the angular difference separating 9 from the specular 
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direction and it) the roughness. Considering two points separated by a horizontal distance 
Ax and a vertical distance Az, we can define the retardation due to these two effects: 

Ac ans = (sin# — sin Q ) Ax , Ac roug h = (cos 9 + cos 9 )Az . 

This allows us to define two characteristic (horizontal) lengths 5 ang and S roug h of the scatter- 
ing system corresponding to the distances between two points of the surface such that Ac ang 
and Ac roug h are equal to the wavelength A. Taking into account the self-afnne character of 
the surface, we get: 

5ang = sin 9- sin 9 ' 5rough = (cos 9 + cos 6 Q fl H * (A) ^ ' 

The coherence length on the surface depends on the relative magnitude of these two charac- 
teristic lengths. For scattering angles close to the specular direction, we have S roug h <C S ang 
and for large scattering angles 5 ang < 5 roug h and the diffuse tail is controlled by the angular 
distance to the specular direction. In general we can evaluate the competition of these two 
effects and their consequences on the scattering cross-section by the simple ratio of the two 
characteristic lengths: 

_ trough _ sin 9 - sin 9 (\\-i/h 

X ~ Sang " (COS 9 + COS 9§)^/^ ' 

We can then describe our scattering system with this unique variable x which takes into 
account the incidence and scattering directions, the roughness parameters of the surface 
and the wavelength. A direct application is the determination of the angular width w of 
the specular peak. The transition between the specular peak and the diffuse tail is simply 
defined by x = 1 which leads to: 

w ~ [2s(A)] 1/if (cos6'o) 1/ ' ff_1 

which is identical to the exact result (21) apart from a numerical constant. Assuming 
that most of the intensity is scattered within the specular peak, we obtain via the energy 
conservation 
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dR s 



~ 1 ~ [2s(A)]~ 1/if (cos 6> ) 1_1/H 



\ <9# 

Forgetting the numerical constants, we can thus rewrite the scattering cross-section as 

\ 86 I S (\y/ H {X) 

When approaching the specular direction we note that S ang diverges whereas 5 roU gh saturates 
at a finite value independent on the angular direction. In this specular direction, the scat- 
tering process is thus controlled by only the latter length and does not depend on the ratio 

X = trough/ Sang- This imposes: 

The argument x being inversely proportional to the quantity s(A) 1 ^ which is nothing but 
a roughness amplitude parameter, the behavior of \& for large arguments can be found 
by matching our expresion with the limit of very smooth surfaces. In this limit a simple 
perturbation approach leads to: 

-2tt 



%)«» 



■(sin 9 — sin 9 ) 
A 



where V is the power density function of the height profile. In the case of a self-afnne profile 
of roughness exponent H, we have V(k) oc k~ l ~ 2H . One can check that this can only be 
consistent with the same power law behavior for 



« X ^ 2H , (x » i). 

V. NUMERICAL SIMULATION RESULTS AND DISCUSSION 

The results obtained in the previous sections were all based on the Kirchhoff approxi- 
mation, and will therefore only be accurate in cases where single scattering is dominating. 
In this section, however, we will therefore no longer restrict ourselves to single scattering, 
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but instead include any higher order scattering process. This is accomplished by a rigorous 
numerical simulations approach which will be described below. This approach will also serve 
as an independent check of the correctness of the analytic results (17), and the results that 
can be derived thereof. Furthermore, it will also provide valuable insight into which part of 
parameter space is dominated by single scattering processes, and thus where formulae (17) 
can be used with confidence. 

The rigorous numerical simulation calculations for the mean differential reflection co- 
efficient were performed for a plane incident s-polarized wave scattered from a perfectly 
conducting rough self-affine surface. Such simulations were done by the now quite standard 
extinction theorem technique [36]. This technique amounts to using Green's second integral 
identity to write down the following inhomogeneous Fredholm equation of the second kind 
for the source function Af(x\u) (see Refs. [39,40]): 

M{x\uo) = 2M Q {x\uo)-2vJ dx' d n G (x, z\x' , z')\ z , =axl) M(x'\Lu). (25a) 

In this equation 

Af(x\u) = d n $(x, z\u)\ z=ax) , (25b) 

where d n = —('(x)d x + d z is the (unnormalized) normal derivative of the total electric field 
$ = E y evaluated on the randomly rough self-affine surface, Af (x\uj) has been defined earlier 
as the normal derivative of the incident field, and V is used to denote the principle part of 
the integral. Moreover, G (x, z\x', z') is the (two-dimensional) free space Greens function 
defined by 

G (x, z\x', z') = itHP |r - r'|) , (25c) 

where r = (x,z), r' = (x',z') and Hq 1 \x) denotes the th order Hankel function of the 
first kind [41]. By taking advantage of Eq. (13a) which relates the scattering amplitude to 
the normal derivative of the total field on the random surface, the scattering amplitude can 
easily be calculated, and from there the mean differential reflection coefficient. It should be 
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noticed that the Kirchhoff approximation used in the previous section to obtain the analyt- 
ical results (17), is obtained from by Eq. (25a) by neglecting the last (integral) term that 
represents multiple scattering. By using a numerical quadrature scheme [42], the integral 
equation, Eq. (25a), can be solved for any given realization of the surface profile ((x). From 
the knowledge of M(x\oo) one might then easily calculate the mean DRC 

Randomly rough Gaussian self-afline surfaces of given Hurst exponent were generated 
by the Fourier filtering method [43] (see Eq. (5)), i.e. in Fourier space to filter complex 
Gaussian random uncorrelated numbers by a decaying power-law filter of exponent —H — l/2 
and thereafter transforming this sequence into real space. The topothesis (or slope) of the 
surfaces were then adjusted to the desired values, £, by taking advantage of Eq. (3). This is 
done by first calculating the topothesy, £ , of the original surface over its total length and 
thereafter rescaling the profile by (£ /£Y~ H , where £ is the desired topothesy. In order to 
having enough statistical information to be able to calculate a well-defined topothesy £o, we 
in fact used a window size slightly smaller than the total length of the surface. 




6 [deg] 6 [deg] 



FIG. 4. A comparison plotted in (a) linear- and (b) linear-log scale, of the mean differential 
reflection coefficient (dR s /d6) vs. scattering angle 6 for a perfectly conducting self-affine surface 
obtained by a rigorous numerical simulations approach (solid lines), and therefore including all pos- 
sible multiple-scattering processes, and the single-scattering results obtained from Eqs. (17) (dashed 
lines). The surface parameters were H = 0.7 and I = 10~ 4 A (s(A) = 0.063) with A = 612. 7nm. 
The angles of incidence for the light were 0° and 50° as indicated in the figure. 
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By the methods just described, we have performed rigorous numerical simulations for 
the mean differential reflection coefficient, (dR s /d0), in the case of a s-polarized plane 
incident wave of wavelength A = 2-kc/uj = 612.7 nm that is scattered from a perfectly 
conducting self-affine surface characterized by the Hurst exponent, H, and the topothesy £. 
For all simulation results shown, the length of the surface was L = 100A, and the spatial 
discretization length was Ax ~ A/10. All simulation results presented were averaged over 
N(; = 1000 surface realizations (or more). Furthermore, in order to check the quality of the 
numerical simulations, both reciprocity and unitarity were checked for all simulation results. 
It was found for all cases considered that the reciprocity was satisfied within the noise level 
of the calculations, while the unitarity was fulfilled within an error of a fraction of a percent. 

In Figs. 4 the mean differential reflection coefficients for a surfaces characterized by the 
parameters H = 0.7 and £ = 10~ 4 A (s(A) = 0.063) are presented. The angles of incidence 
of the light were 8 = 0° and 50° as indicated in the figures. The solid lines represents 
the numerical (multiple-scattering) simulation results while the dashed lines are the (single- 
scattering) prediction of Eqs. (17). As can be seen from Fig. 4a the correspondence is quite 
good between the analytic results and those obtained from the numerical simulations. To 
allow a better comparison for large scattering angles we present in Fig. 4b the results of 
Fig. 4a, but now in a linear-log scale. From this figure it is apparent that for the largest 
scattering angles there are some disagreements between the analytic and numerical results. 
The analytic results tend to overestimate the mean DRC in these regions. This discrepancy 
stems from the fact that multiple scattering is not included in the analytical results. Part of 
the light that according to single scattering would have been scattered into large scattering 
angles are now due to multiple scattering processes, scattered back into smaller 
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FIG. 5. The same as Figs. 2 (single scattering results), but now using a rigorous numerical 
simulation approach (see text for details) that incorporate all higher order scattering processes. 




l[k] 

FIG. 6. The specular peak amplitude, {dR s /d6)\ e=e , and its half width at half maximum, 
w(H,£/X,9q) as a function of the topothesy £. The angle of incidence was in both cases 9q = 0°. 
The solid lines are analytical results obtained from Eqs. (20) and (21) respectively, while the circles 
(amplitudes) and the squares (widths) are obtained from the numerical simulations results shown 
in Fig. 5a. 

angles. This results in smaller values for (dR s /d6) for the largest scattering angles. Since 
the unitarity condition, Jl^ 2 (9R s /d6) d6 = 1, is satisfied for a perfectly reflecting surface, 
this large angle reduction of the mean DRC has to be compensated by an increase for other 
scattering angles. In the case of normal incidence (9 = 0°) say, this increase can be seen 
in the region around \9\ ~ 25° where the numerical simulation results are bigger then the 
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corresponding single-scattering analytical results. The same behavior can be observed for 
an angle of incidence of 50°. 

We give in Figs. 5 the numerical simulation results for five different values of the topothesy 
ranging from i = 1(T 6 A (S(X = 0.016) up to 10~ 2 A (s(A) = 0.25). These multiple-scattering 
results should be compared to the results of Figs. 2 which show the corresponding curves 
obtained from Eqs. (17). The roughness exponent used in the simulations leading to the 
results of Figs. 5 was in all cases H = 0.7, while for the angles of incidence we used 9 = 0° 
(Fig. 5a) and 9 = 50° (Fig. 5b). The height standard deviation as measured over the whole 
length of the surface, L = 100A, was according to Eq. (3) ranging from a(L) = 0.4A for the 
smallest topothesy up to as large as 6.3A for the largest. The fact that we did not really 
use the total length, L, during the surface generation when adjusting the topothesy, but 
instead a slightly smaller fraction of this length, did not seem to affect the height standard 
deviation to a large degree. In fact it was found numerically that the RMS-heights of the 
generated surfaces were only a few percent lower then the one obtained from using Eq. (3) 
and we will therefore in the following use this equation in estimating the RMS-height of the 
surfaces. According to optical criterion these surface roughness correspond to rather rough 
surfaces. In particular one observes from Fig. 5 that in the case of £ = 10~ 2 A a specular 
peak is hard to define at all in the mean DRC spectra. This is a clear indication of a highly 
rough surface and thus a very severe test of our theory. 

To further compare the analytic results derived earlier with those obtained from the 
numerical simulation approach, we in Fig. 6 have plotted the amplitude of the specular 
peaks (circles) {dR s /d9)\ e=eo , and their width (squares) w(H,£/\,9 ), as obtained from the 
numerical simulation results shown in Figs. 5. The solid lines of this figure are the analytic 
predictions for these quantities as given respectively by Eqs. (20) and (21). As can be seen 
from this figure, the analytic predictions are in excellent agreement with their numerical 
simulation counterparts. In particular this confirms the decaying and increasing power-laws 
in topothesy of exponent 1/H — 1 for these two quantities respectively. 

From Eqs. (17) we observe that if we replot the mean DRC times the inverse of the 
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prefactor of the Levy distribution, vs. its argument, all mean DRC-curves corresponding 
to the same Hurst exponent should (within single-scattering) collapse onto one and the 
same master curve. This master curve should be the Levy distribution, £2h(x), of order 
2H. Notice that this data-collapse should hold true for arbitrary values for the angle of 
incidence and topothesy. The failure of such a data-collapse (onto C 2 h) indicates essential 
contributions from multiple-scattering effects. The range of scattering angles where such 
processes are important can therefore be read off from such a plot. Furthermore, since 
the tails of the Levy distribution C 2 h(x) drops off like x" 211 ^ 1 (cf. Eq. (23)) such rescaled 
mean DRC plots can be used to measure the Hurst exponent of the underlying self-afnne 
surfaces for which the light scattering data have been obtained. In order to check these 
predictions for our numerical simulations results, we present in Fig. 7 such a rescaling of the 
data originally presented in Fig. 5b. Only data lying to the left of the specular peak have 
been included, i.e. only data for scattering angles 9 < 9 = 50°. As can be seen from this 
figure the various scattering curves nicely fall onto the master-curve (solid line) in regions 
where single-scattering is dominating. When 
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FIG. 7. Rescaled version of the rigorous numerical simulation results shown in Fig. 5b. Only 
the data corresponding to 6 < 50° are included. In the rescaled coordinates all data (symbols) 
should within the single scattering approximation collapse onto a Levy distribution of order 2H 
(solid line). 



20 



multiple scattering processes start giving a considerable contribution the scattering curves 
start to deviate from this master-curve. This observation could be used in practical ap- 
plications to determine for what regions the scattering is dominated by single scattering 
processes. For the lowest topothesy considered here, £ = 10~ 6 A, a power-law extends nicely 
over large regions of scattering angels — a signature of the diffuse scattering from self- 
affine surfaces. According to Eq. (24) the exponent of this power-law should be —1 — 2H. 
A regression fit to the scattering curve corresponding to the topothesy £ = 10~ 6 A gives 
H = 0.73 ± 0.02, where the error indicated is a pure regression error. The real error is of 
course larger. With the knowledge of the Hurst exponent obtained from the decay of the 
diffuse tail of the mean DRC, we might now, based on the amplitude of the specular peak, 
obtain an estimate for the topothesy of the surface. From the numerical simulation result 
we have that (dR s /d9)\ d=do ~ 17.9 which together with Eq. (20) gives £ = 0.97 • 10~ 6 A, 
where we have used the value found above for the Hurst exponent. These two results fit 
quite nicely with the values H = 0.7 and £ = 10~ 6 A used in the numerical generation of the 
underlying self-afnne surfaces. 

It should be noticed that for the numerical results presented in this paper, we have 
not considered topothesies smaller then £ = 10~ 6 A. However, since lowering the topothesy 
will, as also indicated by our numerical results, favor single-scattering processes over those 
obtained from multiple scattering, the analytic results (17) will trivially be valid for low 
values of the topothesy. This has also been checked explicitly by numerical simulations 
(results not shown). 

So far in this paper we have assumed that the metal was a perfect conductor. Obviously 
this is an idealization, and even the best conductors known today are not perfect conductors 
at optical wavelengths. By relaxing the assumption of the metal being perfectly conducting 
to instead being a good conductor, i.e. a real metal, we are no longer in position to obtain a 
closed form solution of the scattering problem, the reason being that the boundary conditions 
are no longer local quantities. In this latter case we therefore have to resort to numerical 
calculations in any case. In order to see how well our analytic results (17) describe the 
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scattering from real metals (in contrast to perfect conductors) we in Fig. 8 give the mean 
DRC, as obtained from numerical simulations [36], for a self-affme silver surface of Hurst 
exponent H = 0.7 and topothesy £ = 10~ 4 A. We recall that this choice for the topothesy 
corresponds to a rather rough surface where the RMS-height measured over the whole length 
of the surface is a(L) ~ 1.45 A. Furthermore, the angles of incidence were 9 = 0° and 50° 
and the wavelength of the incident light was A = 612.7 nm. At this wavelength the dielectric 
constant of silver is e(uo) = —17.2 + 0.50i [44]. The long dashed lines of Fig. 8 represent the 
predictions from Eqs. (17), and as can be seen from this figure, the correspondence is rather 
good. It is interesting to see that the agreement between the analytical and numerical results 
is of the same quality as the one found for the perfect conductor (see Fig. 4b). This indicates 
that the analytic results given by Eqs. (17) are rather robust and tend to also describe well the 
scattering from a good, but not necessary perfect, reflector. Simulations equivalent to those 
reported for silver have also been performed for aluminum (results not shown) which has a 
dielectric function that is more then three times higher at the wavelength (A = 612. 7nm) 
used here. The conclusions found above for silver also hold true for aluminum. We find it 
interesting to note that such self-affme aluminum surfaces were recently reported to be seen 
for cold rolled aluminum [17]. The Hurst exponents were measured to be H = 0.93 ± 0.03 
and H = 0.50 ± 0.05 for the transverse 
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FIG. 8. The same as Fig. 4b, but now using a real metal (silver) instead of a perfect conductor. 
The value of the dielectric constant of silver at the wavelength of the incident light (A = 612. 7nm) 
was e(uj) = -17.2 + 0.50i. 
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and longitudinal direction, respectively. Before closing this section it ought to be mentioned 
that for real metals the numerical simulations approach based on Eq. (25a), and used above, 
can no longer be used directly. Instead a coupled set of inhomogeneous Fredholm integral 
equations of the second type have to be solved for the electric field, which is non-zero on 
the surface of a real metal, and its normal derivative divided by the dielectric constant of 
the metal. Details about this approach can be found in e.g. Ref. [36]. 

VI. CONCLUSIONS 

We have considered the scattering of s-polarized plane incident electromagnetic waves 
from randomly rough self-afline metal surfaces characterized by the roughness exponent, H, 
and the topothesy, £ (or slope s(A)). By considering perfect conductors, we derive within 
the Kirchhoff approximation a closed form solution for the mean differential reflection coef- 
ficient in terms of the parameters characterizing the rough surface — the Hurst exponent 
and the topothesy (or slope) — and the wavelength and the angle of incidence of the inci- 
dent light. These analytic predictions (written from a Levy distribution of index 2H) were 
compared against results obtained from extensive, rigorous numerical simulations based on 
the extinction theorem. An excellent agreement was found over large regions of parameter 
space. Finally the analytic results, valid for perfect conductors, were compared to numeri- 
cal simulation results for a (non-perfectly conducting) aluminum self-afline surface. It was 
demonstrated that also in this case the analytic predictions gave quite satisfactory results 
even though strictly speaking they were outside their region of their validity. 
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